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Abstract 

We study numerically the directed motion of vortices (antivortices) under an applied ac bias in two-dimensional 
Josephson junction arrays (JJA) with an asymmetrically modulated periodic vortex pinning potential. We find 
that the ratchet effect in large 2D JJA can be obtained using the RSJ model for the overdamped vortex dynamics. 
The rectification effect shows a strong dependence on vortex density as well as an inversion of the vortex flow 
direction with the ac amplitude, for a wide range of high magnetic field around / = 1/2 ( / being the vortex 
density). Our results are in good agreement with very recent experiments by D.E. Shalom and H. Pastoriza [Phys. 
Rev. Lett. 94, 177001, (2005)]. 

Key words: Josephson junction arrays, ratchet effect, periodic pinning 
PACS: 74.25. Qt 



1. INTRODUCTION 

The study of ratchet effects in superconductors 
has attracted great interest in recent years due 
to the possibility of controlling vortex motion [1] . 
However, only very recently nano-micro engineer- 
ing has allowed to successfully fabricate different 
types of asymmetric pinning potentials for the mo- 
tion of flux quanta. Voltage rectification was spec- 
tacularly observed in Nb films with a triangular 
dots array of pinning potential [2] and in Al films 
with square arrays of hole pinning sites [3] . 

Ratchet effect has also been analyzed theoret- 
ically and experimentally in Josephson Junctions 
arrays (JJA), both in the classical and quantum 
regime, particularly for the cases of ID parallel 
JJA, SQUID's or circular arrays [4]. Less atten- 
tion has been paid however to larger two dimen- 
sional JJA, where collective effects are expected to 



play an important role. Indeed, this collective be- 
havior was very recently observed in experiments 
with large square JJA [5]. In this case, improved e- 
lithography techniques were used to modulate the 
gap between the superconducting islands and thus 
an asymmetric and periodic sawtooth potential for 
the vortex-antivortex motion was generated. Volt- 
age rectification was then clearly observed under 
an applied ac current (rocking ratchet), and ana- 
lyzed as a function of the vortex density. In this 
paper we show that the main experimental results 
reported in [5] can be successfully reproduced and 
analyzed in detail by numerical simulations of the 
overdamped RSJ model for a large 2D JJA with 
asymmetrically modulated critical currents. 
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2. MODEL 

We study the dynamics of asymmetrically modu- 
lated 2D JJA using the resistively shunted junction 
(RSJ) model for the SNS (superconductor-normal- 
superconductor) junctions in square two dimen- 
sional arrays [6]. The current flowing in the junc- 
tion between two superconducting islands is the 
sum of the Josephson supercurrent and the normal 
current (see details of the model in [7] ) : 

^W^osin^N+ ^y + ,,(n,t)(l) 

where io is the critical current of the junction 
between the sites n and n + fj, in a square lattice 
[n = (n x ,n y ), = x, y], Rn is the normal state 
resistance and # M (n) = 9(n + ^i) — 9(n) — A^(n) = 
A M 6>(n) — A^(n) is the gauge invariant phase 

difference with A„(n) = §^ J^T^ A ' dl bcin S 
<f>o the flux quantum. The Langevin noise term 
r]^ models the contact with a thermal bath at 
temperature T and satisfies (r/ M (n, t)r}^ (n', t')) = 

^Wn,n'<*(t-0. 

The asymmetry in the underlying periodic pin- 
ning potential is introduced by modulating the 
critical currents Io, and therefore the Josephson 
coupling energy between islands Ej = & I /2tt. 
Experimentally, the ratchet potential is generated 
by modulating the gap between the superconduct- 
ing islands, d, as a sawtooth with a given period, 
p [5]. For SNS junctions, the coupling energy de- 
crease exponentially with d: Ej « exp[— d/£jv], £n 
being the coherence length in the normal metal. 
In our model we make a simplification in order to 
modulate the energy couplings through the criti- 
cal currents. We consider a sawtooth modulation 
of Io where Io = f(n x ) which increases linearly 
from Io min to Io max with each period p. 

An external magnetic field H perpendicular to 
the sample is applied, such that A^ x ^(n) = 
2nf, where / = Ha 2 /&o is the vortex density and 
a is the array lattice spacing. In addition, we apply 
an external ac driving, I = I ac sin(27rwi) in the y 
direction. We take periodic boundary conditions in 
both directions. 

The dynamical equations for the superconduct- 
ing phases are obtained after considering conserva- 



tion of the current in each node (see details in [7]). 
The resulting set of Langevin equations are solved 
using a second-order Runge-Kutta algorithm with 
time step At = O.It,/ (t,/ = 2-kcRmIo/^o) an d 
time integration ranges 2.10 6 tj — 1.10 7 tj after 
the transient. Under ac external currents applied 
perpendicularly to the sawtooth modulation, the 
current-voltages characteristics are calculated as 
well as the vortex structure. The normalization 
used is the following: currents by Io, voltages by 
RnIo, temperature by Io^o/^ks and frequencies 
by {rj)-\ 

We study the vortex rectification on ratchet- like 
square [L x L) JJA, with period p — 8 and linear 
modulated critical currents between Io min = 0.5/o 
and Io max — l-5/o- A wide range of magnetic field 
was explored, from the one corresponding to only 
one vortex in the system (/ = 1/L 2 ) to very high 
vortex densities. We use system sizes from L = 32 
to L = 128. 



3. RESULTS 

It is well known that the velocity of an over- 
damped particle under an applied ac force can 
be rectified by introducing a periodic potential 
with broken reflection symmetry. We can thus ask 
whether this effect, known as rocking ratchet, can 
be still observed for the motion of interacting vor- 
tices in a modulated 2D JJA. Since vortex velocity 
is proportional to voltage, this kind of ratchet ef- 
fect would be observable as a voltage rectification. 
In Fig.l we show clear evidence of the latter ef- 
fect, i.e a net directional vortex motion under the 
applied ac-drive of zero mean. We see that mean 
dc voltage, (Vd c ), as a function of I ac presents 
a maximum for an optimal value of I ac and de- 
creases slowly with increasing I ac . This is the 
qualitative behavior expected for the one particle 
rocking ratchet. To study the role of the vortex 
interactions in Fig. 1(a) examples for different vor- 
tex densities are presented, / = 1/8,1/16,1/32. 
A decrease in the maximum rectification is seen 
while the vortex density increases as well as a shift 
to higher optimal ac amplitude. An important 
feature is worthy of being mentioned here (see in- 
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Fig. 1. Rectified vortex motion: (a) I ac — V characteris- 
tics for low vortex densities. Inset: Comparison with one 
single vortex, (b) I ac — V characteristics for higher vortex 
densities. Simulations performed in 32 X 32 system size, at 
T=0.05 and uj = 1CT 4 . 
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Fig. 2. Rectified dc voltage response dependences: (a) On 
frequency, w, for / = 1/8 and I ac = 0.25. (b) On temper- 
ature, T, for Iac = 0.1 and I a c = 0.25, both examples at 

w = io- 4 . 



set in Fig. 1(a)): the voltage rectification is highly- 
negligible when only one single vortex is present 
in the system, corresponding to a vortex density 
/ = 1/(32 x 32) in this case. A more intriguing 
feature appears when the vortex density is further 
increased. Examples are shown in Fig. 1 (b) . For in- 
stance, an inversion of the vortex flow direction is 
observed. These observations provide evidence for 
the relevant role of vortex interactions. Also notice 
the expected reflection symmetry around f — 1/2 
due to the symmetries with the magnetic field of 
the Hamiltonian representing J J A [8] . From / = 
to / = 1/2 the dissipation arises from the motion 
of "positive" vortices. The opposite occurs for / 
between / = 1/2 to / = 1 where a flux of negative 
vortices (antivortices) gives a voltage response of 
opposite sign. In addition to reproducing the ex- 
periments of Ref. [5], we have observed that the 
instantaneous vortex lattice dynamics behind this 
phenomenon is very complex [9] compared with 
the vortex motion observed in Ref. [2] . 

It is important to accurately determine the 
range of parameters where vortex motion rectifica- 
tion is observable. In agreement with experiment, 
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Fig. 3. Size analysis: (a) Current-voltage characteristics for 
a high vortex density (/ = 3/4) and different system sizes 
L = 32, 64, 128. (b) Maximum and minimum optimum 
mean voltages vs system size for / = 3/4 and / = 1/8. 



3 



0.8 {&) 




0.0 0.2 0.4 0.6 0.8 

I /I 

ac' o 



0.8 
0.6 



1 2 3 4 5 

I [mAl 
ac L 1 

Fig. 4. Rectified voltage vs / and ac drive amplitude. Red 
dotted lines represent (V)dc < and green solid lines, 
(V)dc > 0. (a) Simulations at T = 0.05 and oj = 10~ 4 
in 32 X 32 arrays with p=8. (b) Experiments from [5] at 
T=3.8K and 2-klu = 1kHz. 

we obtain almost zero voltage rectification for all 
ac amplitudes if the frequencies are larger that 
O.OOl(rj)^ 1 due to the fast changes in potential 
slope not allowing vortices to have time to explore 
the asymmetric underlying potential during a cy- 
cle. A plateau of saturation is obtained for very 
low frequencies ( see Fig. 2(a) where an example 
is shown for / = 1/8 and I ac = 0.25). Most of 
the results shown here were also obtained at zero 
temperature as expected for rocking ratchets. It 
is also useful to explore the range of temperatures 
where rectification occurs. An analysis is shown in 
Fig. 2(b) for two values of the ac amplitude which 
are representative of the behavior of almost all the 
vortex densities studied (see also Fig.l): I ac = 0.1 
and I ac = 0.2. 

In Fig. 3(a) we show that our main results at 
high vortex density (/ = 3/4) are independent of 
system size. The key features, such as vortex recti- 
fication and current reversal effects do not depend 
strongly on system sizes. An analysis for the maxi- 
mum and minimum values of the rectification vs L 
is shown in Fig. 3(b) for two vortex densities, / = 
1/8 and / = 3/4 (for more detailed explanation 



on the dynamics at intermediate ac amplitudes see 
[9])- 

Finally, in Fig. 4, contour plots show the occur- 
rence of voltage rectification in the I ac — / plane. 
The plot was generated from (V)dc vs Iac curves 
obtained for different magnetic fields as those 
shown in Fig. 1 and Fig. 3(b). A qualitatively 
good agreement with experiments [5] is obtained. 
A more detailed inspection shows some differences 
comparing both numerical-experimental data. 
The denser and softer experimental contours are 
obtained by continuously changing the magnetic 
field while the simulations were performed for a 
discrete number of vortex densities. 

In conclusion, we find rocking ratchet effects 
in asymmetric modulated 2D JJA using the RSJ 
model for overdamped vortex dynamics. A strong 
dependence on vortex density as well as an inver- 
sion of the vortex flow is found, in good agreement 
with experiments. 
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